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Abstract 



PLh 



We present a measurement of the H i column density distribution function [/(A'^hi)] at 

the present epoch for column densities > 10^'^ cm~^. These high column densities compare 

to those measured in damped Lya lines seen in absorption against background quasars. 

^S) ■ Although observationally rare, it appears that the bulk of the neutral gas in the Universe 

^P i is associated with these damped Lya systems. In order to obtain a good anchor point at 

z = we determine /(A'hi) in the local Universe by using 21cm synthesis observations of 

a complete sample of spiral galaxies. We show that /{Nm) for damped Lya systems has 

O,. changed significantly from high z to the present and that change is greatest for the highest 

' I column densities. The measurements indicate that low surface brightness galaxies make a 

^ • minor contribution to the cross section for H I, especially for A'^hi > lO'^^ cm~^. 
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5^ 1 Introduction 

High column density absorbers seen in the spectra of background QSOs are referred to as Damped 
Lya (DLa) systems if the observed H I column density exceeds the value of A^hi = 2 x lO^*' cm~^. 
The DLa absorption lines are on the square-root part of the curve of growth where damping 
wings dominate the profile and column densities can be determined accurately by fitting the 
line profiles. Wolfe (1995) argues that these systems at high redshifts are gas-rich disks in the 
process of contracting to present-day spiral galaxies. This idea is supported by the fact that the 
characteristic velocity profiles of metal lines and Lyman series lines in DLa systems are similar 
to those of sightlines through spiral galaxies ai z = 0. More recently, detailed modeling of DLa 
absorption profiles by Prochaska & Wolfe (1998) has shown that the DLa systems are consistent 
with rapidly rotating, thick disks. Note however that alternative models, like protogalactic clumps 
coalescing into dark matter halos (Haehnelt et al. 1997, Khersonsky & Turnshek 1996), can also 
explain the kinematics. The cosmological mass density of neutral gas in DLa systems at high 
redshift is comparable to the mass density of luminous matter in galaxies at z = (e.g. Lanzetta 
et al. 1995). 

One of the best known statistical results of the study of QSO absorption line systems is the 
column density distribution function (CDDF) of neutral hydrogen. The function describes the 



chance of finding an absorber of a certain H I column density along a random line of sight per unit 
distance. An observational fact from high-z Lya studies is that the differential CDDF [/(A^hi)] 
can be described by a single power law of the form /(A'hi) oc Ahi", where a ^ — 1.5 over ten 
orders of magnitude in column density (e.g. Tytler 1987, Hu et al. 1995) from 10^^ cm~^ (Lya 
forest) to 10^2 cm'^ (DLa). 

An integration over the distribution function gives the total cosmological neutral gas density 
as a function of redshift. The H I gas density relates to /(Ahi) as ^hi oc Jj^^ N}iif{Nui)dNui and 
it is readily seen that i^Hi(AHi) oc A'^'^ if a = —1.5 and N2 ^ A^i. This implies that although 
the high column density systems are observationally rare, they contain the bulk of the neutral 
gas mass in the Universe. Because so few DLa systems are known {^ 80), the uncertainties on 
JIhi and the CDDF for high column densities are large, especially if the measurements are split 
up into different redshift bins. But following the CDDF as a function of redshift is certainly 
very important in constraining models of complicated physical processes like star formation or gas 
feedback to the interstellar medium. 

There are several reasons why the determination of /( Ahi) at the present epoch is difficult. Due 
to the expansion of the Universe the expected number of absorbers along a line of sight decreases 
with decreasing redshift, the Lya line is not observable from the ground for redshifts smaller than 
1.65, and starlight and dust in the absorbing foreground galaxies hinder the identification of the 
background quasars. Gravitational lensing may also play a role as it can bring faint quasars into 
the sample which otherwise would not have been selected (e.g. Smette et al. 1997). 

At the present epoch the largest repositories of neutral gas are clearly galaxies. No instance 
of a free-fioating H I cloud not confined to the gravitational potential of a galaxy has yet been 
identified. It is therefore justified to use our knowledge of the local galaxy population to estimate 
the shape and normalization of the CDDF. 

2 How to determine /(A^m) at ^ = 0? 

A simple but illustrative and instructive method is to take the analytical approach. This is 
illustrated in Figure |I[ Here we represent the radial distribution of the neutral hydrogen gas in 
galaxies by both an exponential and a Gaussian model. The differential cross sectional area of an 
inclined ring with a column density in the range A^ to N+dN is given by dIl{N, i) = 2'Kr{N)dr cos z, 
where r{N) is the radius at which a column density A^ is seen, and i is the inclination of the ring. 
We assume that the luminosity function 0(M) of the local galaxy population can be described 
by a Schechter function as indicated in the upper right panel of figure 1. The local f{N) can be 
derived from (p{M) and the area function c?S(A) by taking the integral 

Hq dN 

where the subscript i indicates an average over all inclinations. To evaluate this integral, the area 
function, or more generally the radial H I distribution, needs to be related to M. Here we adopt 
the relation log Mhi = A + BMb (following Rao & Briggs 1993) and assume that the central gas 
surface density in disks is not dependent on morphological type or luminosity. The resulting f{N) 
for both models is shown in the lower left panel. The integral H I gas density in /iioogcm"^ as 
a function of column density is shown in the lower right panel. This function can be calculated 
with Phi(A) = miiN-^f(N)dN, where mu is the mass of the hydrogen atom. 

The Gaussian models yield a CDDF of the form /(A) oc A^", where a = —1 for A^ smaller 
than the maximum column density seen in a face-on disk (Amax) and a = —3 for the higher values 
of A^. The exponential model gives a smoother function. The logarithmic slope is approximately 
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Figure 1: Illustration of simple models for the CDDF. Upper left: Gaussian and exponential 
models for the radial distribution of neutral gas density in galactic disks. Upper right: Schechter 
function describing the local luminosity function. Lower left: Resulting CDDFs for the Gaussian 
and exponential radial profiles. Lower right: Integral neutral gas density in the local Universe as 
a function of column density for both models. The effect of disk truncation is indicated by the 
dashed line. 



— 1.2 around A^ = 10^° cm~^, slowly changing to —3 at higher column densities. In fact, it was 
shown already by Milgrom (1988) that f{N) oc N~'^ for A^ > A^max for any radial surface density 
distribution. The lower right panel clearly illustrates that an overwhelming part of the total H I 
mass in the local Universe is associated with column densities close to A^max- 

In addition to these simple models we also show the effect of disk truncation on the CDDF. The 
thin dashed line illustrates a Gaussian disk truncated at A^hi = 10^^'^ cm~^, the level below which 
photo-ionization by the extragalactic UV-background is normally assumed to be important (e.g. 
Corbelli & Salpeter 1993, Maloney 1993). It appears that this truncation only seriously affects 
the CDDF below A^hi = 10^^'^ cm~^. No significant changes occur at higher column densities. 

A more reliable method than this analytical approach is to determine /(A^hi) by using observed 
H I distributions. 21cm maps of nearby galaxies routinely reach sensitivity limits comparable to 
column densities that typify DLa absorbers. It therefore seems natural to calculate /(A'hi)(2; = 0), 
simply by adding cross sectional areas as a function of A^hi for a large sample of galaxies for which 
21cm synthesis observations are available. However, this approach is complicated by the fact 
that there is an enormous variation in sensitivity and angular resolution of the 21cm maps, and 
the problem of choosing a fair and complete sample of galaxies. Most galaxies that have been 



studied extensively in the 21cni line were selected on having either a large H I diameter, so that 
the rotation curve can be sampled out to large galactocentric radii, or on having peculiarities 
such as polar rings or warps. Thus, most samples for which 21cm synthesis data exist are not 
representative of the galaxy population of the local Universe and would likely be biased against 
dwarf and low surface brightness galaxies. 

3 The Ursa Major Cluster 

The Ursa Major cluster of galaxies, studied extensively by Verheijen (1997), forms an ideal sample 
for an unbiased study of the column density distribution function. The Ursa Major cluster, 
at a distance of 15 Mpc, is different in many respects from famous, equally distant clusters 
like Virgo and Fornax. Ursa Major's member galaxies show no concentration towards a central 
condensation and their velocity dispersion is exceptionally low, approximately 150 km/s. The 
estimated crossing time is half a Hubble time and hence the galaxies do not seem to be seriously 
affected by tidal interactions. In addition to this, there is a predominance of late type galaxies and 
the morphological mix of galaxies is indistinguishable from that in the field. This combination 
of properties implies that the Ursa Major cluster is in fact an overdensity of galaxies and not 
comparable to classical clusters of galaxies. This justifies the use of the Ursa Major cluster for 
the study of the shape of the CDDF of neutral hydrogen in the local Universe. 

The Ursa Major cluster as defined by Tully et al. (1996) comprises a volume of 80 Mpc^, 
within which 80 galaxies are identified to date. For a complete sample of 62 galaxies intrinsically 
brighter than the Small Magellanic Cloud {Mb = —16.5™) 21cm synthesis observations have been 
performed with the WSRTQ. H I has been detected by the WSRT in 49 galaxies. Details on 
observations and data reduction are described in Verheijen (1997). 

An obvious advantage of using the UMa sample for this study is that all the member galaxies are 
at approximately the same distance. Therefore, the spatial resolution of the synthesis observations 
are constant for the whole sample. This simplifies the problem of assessing the influence of 
resolution on the determination of the CDDF and the comparison with the CDDF at high redshift. 

The shape of the column density distribution function is determined by counting in each H I 
map the number of pixels per logarithmic bin of 0.1 dex in column density. The solid angle covered 
by pixels of a certain column density is then determined by multiplying the number of pixels with 
the angular pixel size which varies slightly from galaxy to galaxy. 

The disadvantage of using a galaxy sample taken from a clear cosmic overdensity is that the 
CDDF is not automatically normalized. If we would naively assume that the Ursa Major cluster 
is a representative part of the nearby Universe, we would overestimate the normalization of the 
CDDF by roughly a factor of 12. This factor is obtained by comparing the H I mass function of 
the cluster with that of the field galaxy population (Zwaan et al. 1997). The shape of the Ursa 
Major mass function is indistinguishable from that of the field, but the normalization, 9*, is larger 
by a factor of ~ 12. Ideally, one would use a sample of galaxies with well understood selection 
criteria so that the normalization would occur automatically. Unfortunately, there are no such 
samples available for which H I synthesis observations with sufficient angular resolution have been 
performed. The HIPASS survey, a blind 21cm survey of the whole southern sky, will eventually 
yield a suitable galaxy sample for this purpose, if a representative subsample is followed up with 
the ATCA to obtain high spatial resolution maps. 

There are several methods for normalizing the UMa CDDF. By assuming a local luminosity 
function (LF) or H I mass function (HIMF), each galaxy could be given a weight according to its 

^Thc WSRT is operated by the Netherlands Foundation for Research in Astronomy (NFRA/ASTRON), with 
financial support by the Netherlands Organization for Scientific Research (NWO) 



absolute magnitude or H I mass. However, this method introduces extra uncertainty in the derived 
CDDF, due to uncertainties in the exact shape and normahzation of the LF and the HIMF. Our 
preferred method of normahzing the CDDF is to scale the complete function, not the individual 
contributors to it. This can be achieved by scaling the integral H I mass density that is contained 
under the CDDF: 

r-A^max Ho 



PHI = / mnN^f{N)dN. (2) 

By means of a blind 21cm survey Zwaan et al. (1997) determined pm = 5.8 x lO^/iioo MoMpc"^, 
a result that is in excellent agreement with earlier estimates based on optically selected galaxies. 
Note that dependencies on Hq disappear in the final specification of the CDDF. 

4 The Column Density Distribution Function 

Figure |] shows the CDDF determined from the 21cm observations of the Ursa Major sample. 
From left to right the function is shown for three different resolutions of the H I maps: 15", 30", 
and 60". The solid line is the determined CDDF; the dashed lines indicate the quality of the 
measured column densities. Each pixel in the H I maps has an estimate of the signal to noise level 
assigned to it. In the determination of the CDDF we calculated an average S/N level for each bin 
in column density by averaging the S/N ratios for the individual pixels. The dashed lines show 
the average la errors on the column densities and should be interpreted as horizontal errorbars. 
Nonetheless, they clearly overestimate the real uncertainties on the CDDF as many pixels are used 
in each bin (2500 independent beams for full resolution). The lines merely serve as an indicator 
of the quality of the measurements at each resolution. The thin solid line represents the CDDF 
for a Gaussian model, where /(A^m) oc A^hi~^ for A^hi < 10^^ cm~^ and /(A^m) oc A^hi~^ for higher 
column densities. 

When comparing the CDDFs at different resolution, it appears that the highest resolution 
maps yield the smoothest CDDF. This occurs because there the measurements of column density 
have the lowest S/N ratios. The low resolution (but high S/N) CDDF is in excellent agreement 
with the Gaussian model for 10^° < N < 10^^'^, but for higher column densities, the measured 
curve drops below the model since high column density peaks are smeared away. Going to higher 
resolutions leads to better agreement between the measured curve and the model for the highest 
A^Hi; and at 15" resolution the CDDF follows the iVni"^ distribution up to 10^^'^ cm~^. 

Besides beam smearing two other effects can cause a deviation from the Nm'^ function. Firstly, 
the calculations assume that the gaseous disks are infinitely thin. Observations show that the 
bulk of the H I indeed resides in a thin layer with axis ratio < 0.1 (Rupen 1991). The thin 
disk approximation is therefore valid for moderately inclined disks. However, the highest column 
densities in the models arise in highly inclined thin disks. A small degree of puffiness will prevent 
these high column densities from being observed. The second effect is H I self absorption. The 
theoretical calculation of the CDDF is based on the assumption that the optical depth of the 
neutral gas layer is negligible. Column densities much higher than the maximal column density 
in a face-on galaxy can only be seen in a highly inclined disk when the gas is optically thin. It is 
remarkable that the full resolution CDDF follows the f{N) oc A^hi~^ line up to A^hi = 10^^'^, well 
above the value where H I self absorption is normally assumed to set in. For example, Dickey & 
Lockman (1990) calculate that an H I cloud with T = 50K and an FWHM velocity dispersion of 
10 kms~^ becomes optically thick (r = 1) at column densities A^ = 10^^ cm~^. 

Also shown in Figure ^ are the measurements of /(A^hi) at high redshifts as determined by 
Storrie-Lombardi et al. (1997). We choose not to split up their high-z sample in different redshift 
bins in order to get reasonable signal to noise. The median redshift of the total DLa sample is 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
"\ 15" 


1 


1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 

30" 


1 


1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 
60" 


1 


^ -22 

rg 

E 

o 




- 


-" ^ - -"~" 


■-^^5^ i-i— 1 


- 


^-^-^ 




- 


S -24 
o 


- ^^ 


- 


- 


\^ 


- 


- 


Y^ 


- 


-26 


1 1 1 1 1 1 1 1 1 1 1 1 1 I'll 


\" 


1 1 1 


1 1 1 1 1 1 1 1 1 1 li 1 


\" 


1 1 1 


1 1 1 1 1 1 1 1 1 1 ii 1 


\" 



19 20 21 22 19 20 21 22 19 

log N„ [cm-2] 



20 



21 



22 



Figure 2: The column density distribution function at z = 0. From left to right the function 
is shown for three different resolutions of the H I maps: 15", 30", and 60". The thick solid line 
is the CDDF as measured in the Ursa Major cluster. The dashed lines indicate the average 
la uncertainties in the column density determinations. For reference, a Gaussian model with 
^max = 10^^ cm~^ is also shown as a thin solid line. The points indicate the H I CDDF for 
damped Lya systems at high redshift {(z) ~ 2.5) taken from Storrie-Lombardi et al. (1997). Note 
the strong evolution, especially at the highest column densities. 



z = 2.5. A value of go = 0.5 has been used here. Lower values of go would not significantly 
change the slope of /(A^hi) but would decrease the normalization by approximately a factor of 2. 
Strong redshift evolution of the CDDF from z = 2.5 to the present is apparent. The intersection 
cross-section for H I column density < 10^^'^ cm~^ has decreased by a factor of 6 (factor 3 for 
go = 0) from z = 2.5 to 2; = 0. Higher column densities show a larger decrease, the evolution 
accounting for a factor of 10 (5 for go = 0). Lanzetta et al (1995) report still stronger evolution of 
the higher column densities for higher redshift, although the highest column densities suffer from 
small number statistics and the effect is hardly seen by Storrie-Lombardi et al. (1997). The strong 
evolution of the higher column densities can be understood if gas consumption by star formation 
occurs most rapidly in regions of high neutral gas density (Kennicutt et al. 1994). 

Rao & Briggs (1993) evaluated the CDDF at the present epoch by analyzing Arecibo observa- 
tions of a sample of 27 galaxies with optical diameters in excess of 7'. Double- Gaussian fits to the 
observed radial H I distribution were used to calculate /(Ahi). The disadvantage of this method 
is that the Gaussian fits automatically introduce the Ahi~^ for low Ahi and Ahi~^ for high Ahi- 
In the present study no modeling has been applied. The location of the change of the slope and 
the normalization are in excellent agreement between Rao & Briggs' work and the Ursa Major 
determination. 



5 Contribution of Low Surface Brightness Galaxies 

It has been argued in the literature that low surface brightness (LSB) galaxies might contribute 
a considerable H I cross section. In particular, hinder (1998) explores a scenario in which the 



outskirts of galaxies are responsible for most of the cross section for low column density neutral 
gas (A^Hi < 10^°'^ cm~^). She concludes that Lya absorber counts at low redshifts can be explained 
if LSB galaxies of moderate absolute luminosity with extended low density gas disks are included 
in the analysis. Contrary to this view, Chen et al. (1998) claim that extended disks of luminous 
galaxies can account for most of the observed Lya lines below A^hi = lO^'''^ cm~^. The contribution 
of dwarf and LSB galaxies to the cross section for high column density H I is also unclear. For 
instance, Rao & Turnshek (1998) show that there are no luminous spiral galaxies in the vicinity 
of the quasar 01 363 in which spectrum they identify two low- 2; DLa systems. 

Here we evaluate the contribution of LSB galaxies to the cross section for high column density 
gas at 2; = 0. First we have to address the problem of completeness. The Ursa Major sample 
is essentially a magnitude limited sample. Selection effects against LSB galaxies are therefore 
to be expected. Tully & Verheijen (1997) discuss the completeness of the sample by plotting 
the observed central surface brightness against the exponential disk scale length. Theoretical 
approximations of the visibility limits seem to describe the boundaries of the observed sample 
satisfactorily. We apply the same visibility limits to the H I selected galaxy sample of Zwaan et 
al. (1997) to estimate what fraction of the H I mass density in the Ursa Major cluster could be 
missed in the present study. It appears that galaxies below the optical detection limits of our 
Ursa Major sample probably contain 10% of the total H I density of the cluster. Not surprisingly, 
most of this missed H I density resides in LSB galaxies. Following Tully & Verheijen (1997), the 
separation between LSB and HSB galaxies is made at an extrapolated central surface brightness 
of 18.5 magarcsec"^ in the i^'-band, which roughly compares to 22.0 magarcsec"^ in the i?-band. 

Figure ^ illustrates the contribution of LSB galaxies to the cross section for high column density 
gas, relative to the total galaxy population. We have corrected for the incompleteness by adding 
extra cross section for the LSB galaxies, equally over all column densities, in such a way that 
the mass density in these galaxies increases by an amount equal to 10% of the total H I density. 
The left panel shows the CDDF, the right panel shows the cosmological mass density of H I as 
a function of column density. The full resolution data are used. LSB galaxies do not make a 
significant contribution to the cross section for column densities higher than A^hi = 10^^'^ cm~^. 
Below that value they are responsible for approximately 25% of the cross section. The right panel 
shows that LSB galaxies make a minor contribution to the local neutral gas density, a conclusion 
very much in concordance with the results of Briggs (1997), Zwaan et al. (1997), Sprayberry 
(1998) and Cote et al. (1998). 

6 Conclusions and Discussion 

We have used the present knowledge of the nearby galaxy population to estimate the H I column 
density distribution function at 2; = 0. It is shown that /(A^hi) undergoes strong redshift evolution 
from z ~ 2.5 to the present, especially at the high column densities. The observed evolution in 
fiNui) critically depends on whether the census of H I in the local Universe is complete. Surveys 
in H I and the optical indicate that the density of visible light and neutral gas is dominated by 
luminous, high surface brightness galaxies. The H I surveys routinely reach column density limits 
much lower than what is required to detect the z = counterparts of DLct systems. Since H I mass 
functions pubhshed to date typically lose sensitivity below Mhi = lO^M©, the region of parameter 
space still open to hide a large amount of high column density gas is that of low H I masses. 
Observations to measure the space density of these small H I masses (H I clouds and extreme LSB 
dwarf galaxies) and to evaluate to what extent they contribute to the H I density and the CDDF 
of the local Universe are important next steps. 
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Figure 3: The column density distribution function {left) and H I mass density (right) for high and 
low surface brightness galaxies. High surface brightness galaxies dominate the H I cross section, 
especially for column densities > 10^^'^ cm~^. 
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